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In honor of Peter Olver's 60-th birthday 
Abstract 

We show how Noether conservation laws can be obtained from the particle relabelling symmetries in the Euler- 
Poincare theory of ideal fluids with advected quantities. All calculations can be performed without Lagrangian 
variables, by using the Eulerian vector fields that generate the symmetries, and we identify the time-evolution 
equation that these vector fields satisfy. When advected quantities (such as advected scalars or densities) are 
present, there is an additional constraint that the vector fields must leave the advected quantities invariant. We 
show that if this constraint is satisfied initially then it will be satisfied for all times. We then show how to solve 
these constraint equations in various examples to obtain evolution equations from the conservation laws. We also 
discuss some fluid conservation laws in the Euler-Poincare theory that do not arise from Noether symmetries, 
and explain the relationship between the conservation laws obtained here, and the Kelvin-Noether theorem given 
in Section 4 of Holm, Marsden and Ratiu, Adv. in Math., 1998. 
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1 Introduction 

As Noether did in her famous paper |21j . we are dealing with invariant variational principles. This subject has a vast 
literature and has been a favorite topic for Peter Olvcr, to which he returned many times [51 12"2"1 12"31 12"41 12"51 [2^1 |2"U] . 
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A Lie group transformation that leaves the Lagrangian invariant in Hamilton's principle is called a variational 
Lie symmetry. The correspondence between variational Lie symmetries and conservation laws for Euler-Lagrange 
equations is completely determined by the Noether's First Theorem QjJJ [21] [26]. Namely, every variational Lie 
symmetry yields a conservation law0 Our main goal here is to identify explicitly in terms of Eulerian observables 
the vector fields of the relabelling symmetry transformations under the Lie group G of smooth invertible maps that 
are responsible for some of the well-known conservation laws in the Euler-Poincare theory of fluids with advected 
quantities |15j . In particular, we treat a few hands-on examples in fluid dynamics that recover some famous formulas 
such as helicity of fluids, Ertel's potential vorticity in geophysical fluid dynamics (GFD) and Chandrasckhar's cross- 
helicity for magnetohydrodynamics (MHD). We also discuss the relation of the classical Noether's Theorem with the 
Kelvin-Noether circulation theorem from the Euler-Poincare theory of ideal fluids with advected quantities in [T5] . 
In addition, we discuss conservation laws in the Euler-Poincare theory that do not arise from Noether symmetries. 
Finally, we discuss some applications of Noether's Theorem in image registration problems. 

It seems that every theoretical physicist and many mathematicians eventually feel compelled to write a paper 
about Noether's Theorem. Previous influential papers along similar lines about Noether's Theorem in fluid dynam- 
ics related to the directions taken here include [HJ EQl [12] (TJ (27] [28] El [18] and of course references therein. 

The main content of the paper is: 

1. Section [5] briefly summarises the Euler-Poincare formulation of ideal fluid dynamics with advected quantities. 
In particular, we summarise several simple but useful theorems that are available for studying how the Noether 
theorem associates variational Lie symmetries with conservation laws for fluids. 

2. Section[3]uscs these theorems in a sequence of examples that derive several of the most well-known conservation 
laws for ideal fluids in the Euler-Poincare formulation |15j . 

3. Section [4] points out that not all fluid conservation laws follow from Noether's theorem, by considering the 
counterexample of magnetic helicity for MHD. It also makes a connection between Noether's theorem as 
discussed in this paper, and the Kelvin-Noether circulation theorem discussed in |15j . 

4. Section [5] discusses some numerical issues and applications of these ideas outside of fluid dynamics. Section [5] 
also raises topics for future research inspired by Lie symmetries and Noether's theorem. 

2 Formulation 

We begin by laying out the assumptions that underlie the Euler-Poincare formulation. These are the following. 

1 . There is a right representation of the action of a Lie group G on its tangent space TG and on the vector space 
V. The action on TG x V is denoted by concatenation on the right, as (v g , a)h = {v g h, ah) for g,h G G. 

2. The Lagrangian L : TG x V — > R is right G-invariant. 

3. In particular, if ao G V, define the Lagrangian L ao : TG — > R by L ag (v g ) = L(v g ,ao). Then L ao is right- 
invariant under the lift to TG of the right action of G ao on G, where G ao is the isotropy group of ao. 

4. Right G-invariance of the Lagrangian L permits us to define a reduced Lagrangian I : g x V — > R by 



Conversely, this relation defines for any I : g x V — > R a right G-invariant function L : TG x V — > R. 

5. For a curve g(t) G G, let u(t) := g(t)g(t)~ l and define the curve a(t) G V obtained from the action G x V — > V 
as the unique solution of the linear differential equation with time dependent coefficients 



For fluids, the Lie group G = Diff(R 3 ) is the group of diffcomorphisms of three-dimensional space. This is the 
Lie group of smooth invertible maps defined on R 3 and with smooth inverses^ At time t, the curve g(t) defines 

1 Noether's celebrated paper 1211 contains two major theorems. The present paper discusses only the first of these theorems. For 
good discussions of the second Noether theorem, see e.g. 1211 126 4 19 18]. 

2 Strictly speaking, G = Diff (M 3 ) denotes the connected component at the identity of the diffcomorphisms. Its Lie algebra comprises 
the right-invariant vector fields on R 3 , denoted as £(R 3 ). 



l(v g g 1 ,a Q g x ) = L(v g ,a ). 



(1) 




(2) 



with initial condition a(0) = ao and Lie derivative C u u). 
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the mapping from a reference configuration (known as label space) to the physical domain so that x(t) = g(t)xo 
describes Lagrangian particle trajectories for each label xq. 

Definition 1. The solution a(t) = ao<?(t) -1 of equation is called an advected quantity for fluids, and the 
right-invariant vector field u(t) := gg(t) -1 G X(IR 3 ) is called the Eulerian, or spatial, fluid velocity. 

Remark 2. Examples of advected quantities include the extensive thermodynamic properties that are carried by 
fluid elements such as their heat and mass. Equation ^ means physically that along the flow g(t) of the vector 
field u(t) the fluid elements are to be regarded as closed thermodynamic systems that do not exchange heat and 
mass with their neighbours. 

Some particular examples of advected quantities that we discuss in this paper are: 
f . Scalar fields (0-forms) a(t) = s that satisfy: 

(dt + £ ui t))s= (d t + u- V)s = 0. 
In geophysical models scalar advected quantities includes buoyancy due to heat and salinity. 

2. Density fields (volume forms) a(t) = pdV that satisfy 

(d t + L u{t) ) pdV = (d t p + V • (up)) dV = 0. 
This type of advected quantity is used for the fluid density, or layer depth in shallow water models. 

3. Flux fields (2-forms) a(t) = B ■ dS that satisfy 

(dt + C u(t) ) B ■ dS = (d t B - curl(w x B)) ■ d S = 0. 

This type of advected quantity is used, e.g., for the magnetic flux in magnetohydrodynamics. 

For more discussion of advected quantities, see [TS]. The back-to-labels map that specifies the label of the fluid 
parcel currently at a given spatial position would also be an advected quantity. However, in this paper, we shall 
restrict ourselves to dealing only with Eulerian observables, and the particle label is not observable at any given 
Eulerian point in a fluid flow. 



2.1 Euler-Poincare theorem with advected quantities 

Here we review the approach presented in |15j to obtaining the variational equation of motion, known as the 
Euler-Poincare equation, for general reduced Lagrangians l(u, a) with advected quantities. 

Hamilton's principle SS — for S — J l(u,a) dt with the reduced Lagrangian defined in equation {TJ) may be 
expressed either abstractly as 

= 6S = d l(u, a) dt = jf 1 ^ , 5u^ + ^ , Sa^ dt , (3) 

where angle brackets denote appropriate pairings, or cquivalently in coordinates with Lagrangian £(u, a) dV 

= SS = S [ 1 [ £(u,a)dVdt= [ ' [ ( ^ ■ 5u + ^ So] dV dt, (4) 
J t Jt> J t Jv \Su 6a J 

where T> is the spatial domain with boundary dV on which the fluid velocity has no normal component; that is, 
u • n = 0. The expressions G X* and J£ G V* are variational derivatives in u and a, respectively. We ensure 
that variations in g honour the boundary conditions, by defining Sg = w o g, in which w is a vector field whose 
components satisfy w ■ n = on dT>. 

In the remainder of the paper, we will find it convenient to use a hybrid notation that passes freely between the 
abstract notation and the more explicit coordinate notation, as in equations ([3]) and (j4]). We believe this hybrid 
notation, whose meaning will always be clear from the context, will appeal to a wider readership than the abstract 
notation. Conversely, we will sometimes find that the calculations we need to perform are written more directly in 
the abstract notation using the language of differential forms. 
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The infinitesimal transformations for u and a are |15j 

5u = w — ad u w := w + [it, w] , 5a = —C u 



(5) 



Here a denotes any quantity that is advected with the flow, e.g. scalar tracers s, densities pdV etc. The linear 
operator on u>, ad u , is defined in terms of [u, w], which is the commutator (Lie bracket) of the vector fields u and 
w in 3C(R 3 ). Furthermore, we seek the stationary point 5S = in Hamilton's principle above, subject to Sg = at 
the endpoints t = t$ and t = ti; hence, we also require w = Sg g~ l to vanish at the endpoints. 
Substitution in ((3|) now yields 







to Jt> 



d 51 
dt Su 



ad* 



51 51 
5u 5a 



o a 



where ad* is the dual operator to ad u defined by 

v ■ ad* mdV 



v 



wdVdt 



m ■ &d u vdV, 



5l_ 

5u 



wdV 



(0 



v 



for all vector fields v, and whose explicit formula in components is 

ad* m = V • (u ® m) + {Vu) T m . 

This formula also happens to match the components of the Lie derivative for one- form densities |15j . 

C u (m ■ dx (g) dV^) = (V • (u ® to) + (Vw) T m) • dx ® dV, 

with line element dx and volume element dV. Notation for the diamond operation (o) has also been introduced 
in equation (|6]). The diamond operation is defined by 



v 



5l_ 

5a 



o a 



>dV := 



51_ 

5 a 



(-C w a)dV. 



(7) 



Vanishing of the first term in 
equation, 



for variations that are otherwise arbitrary now produces the Euler-Poincare (EP) 



— — + V 

dt 5u 



51 



51 51 



(8) 



The EP equation in ((SJ) is completed as an evolutionary system by including the equation of motion (J2|) for the 
advected quantities, a. 



Noether's Theorem for Euler-Poincare with advected quantities We consider symmetries of the action 
S = J l(u,a) dt that are obtained by infinitesimal transformations of the form 5g = r/ o g for a vector field to 
Consequently, we have the infinitesimal transformations 

5u = r\ + [u, rj] , 5a = —£ v a . (9) 

If the vector field i] generates symmetries of the Lagrangian, then Hamilton's principle SS = implies that 

f fl f ( d 51 „ / 51 \ sT 6l 51 \ r f 51 1T 1 41 

= -/ / — — +V- u<8>— +(Vu) T — - — oa )-ndVdt+ / —-ndV . (10) 
Jt a Jv\otou \ buj 5u 5a J [J v 5u J. 

^ v ' 

= 

Here, the term in the time integral vanishes for solutions of the Euler-Poincare equations, and we are left with the 
endpoint terms for arbitrary to and t\. This implies Noether's theorem. Namely, a conservation law is associated 
with each vector field n that generates a symmetry of the Lagrangian [21j . These considerations prove the following. 

Theorem 3 (Noether theorem for EP). Each symmetry vector field rj of the EP Lagrangian (0) for infinitesimal 
transformations given by {Pj] corresponds to an integral of the EP motion equation @) satisfying 

HZ-"*- 9 - 

for an appropriate inner product. 
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2.2 Relabelling symmetries 

Let us now consider how to derive the vector fields rj for the symmetry transformations in Noether's theorem in the 
case of fluids in the Eulerian representation. These symmetry transformations are called relabelling symmetries. 
They arise from the assumed right invariance of the EP Lagrangian l{u, a) under the group G of diffcomorphisms 
(the Lie group of smooth invertible maps with smooth inverses). The Eulerian velocity u(t) := ggit)^ 1 G X(R 3 ) is 
right-invariant under this action and therefore it does not change under relabelling transformations. This invariance 
implies the following evolution equation for the vector field r\: 

6u = fi + [u,r)} =0, (12) 

where the bracket [ • , • ] denotes commutation of vector fields. If a set of advected quantities {a} exists, then the 
vector fields r\ for the symmetry transformations must also satisfy the additional conditions that 

5a = — C v a = , (13) 

for each advected quantity a. 

When there are no advected quantities present (as in the case of EPDiff [14] , for example) equation (|TT]) simply 
recovers the equation for conservation of momentum, as one sees from the following direct computation: 

d /SI 

= diV^ 

d SI \ I 51 d 
dt Su' J \5u dt 
d SI \ I 51 



1% V] 

at du I \bu 
A 51 \ I 51 , 

dtfc'7 + W u?7 

d 51 „ 51 
— — + ad u — ,7? 
dt du du 



In this computation, the angle brackets ( • , • } denote the L 2 pairing X* x X — > K between the vector fields and their 
L 2 duals, the 1-form densities. 



2.3 Theorems for advected quantities 

We now develop general results for the case where one or more advected quantities are present. This requires 
determining whether all the conditions in © can be satisfied simultaneously. We shall conclude that, if they are 
satisfied initially, then they are satisfied for all times t, due to the commutative properties of Lie derivatives. This 
will enable us to derive conservation laws in various cases in the rest of the paper. 

Theorem 4 (Commutator). For any pair of smooth time- dependent vector fields u(t), n(t) £ X and for any a(t) G V 
the following commutation relation holds among Lie derivatives, 

[d t + C u{t) , C v ( t) ] a(t) = C^ + [ u ^a(t) . (14) 

Proof. For any a(t) G V, one computes by the product rule for Lie derivatives that 

(d t + C u ( t )) C n a(t) = C(^ + [ UjV ])a(t) + C v (d t + £«(*)) a(t) . 

Hence, the commutation relation in (|14p holds, and because a(i) G V is arbitrary, this implies the Lie derivative 
commutation relation 

[d t + £„(t) , C v ] = C(r) + [ Ui7j ]) . (15) 

□ 
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Under the assumption that the variational vector field n satisfies the time-evolution equation (|T^|) required for 
a relabelling symmetry transformation, one finds the following commutator theorem. 

Corollary 5 (Symmetry). If a vector field r\ satisfies equation for an infinitesimal relabelling symmetry, then 
the Lie derivative C v commutes with the evolution operator, (pt +C U U)), 

[dt + A*(t) i A,(t)] a (*) = f° r V + K v] = • ( 16 ) 

Proof. This symmetry corollary follows by inserting equation f|12[) into the commutation relation in equation fj 14|) . 

□ 

Theorem 6 (Ertcl theorem). If the quantity a is advected as in equation {![) and the vector field r] satisfies equation 
for an infinitesimal relabelling symmetry, then L n a is also advected. 

Proof. By equations and (THZ1) one finds the advection relation for £ n a, 

(d t + C u{t) ) C n a(t) = Cr, {d t + C u[t) ) a(t) = , (f 7) 

as a result of the condition ([2]) satisfied by advected quantities. □ 

Consequently, if C v a = in equation (|13|) holds initially, then it continues to hold under the EP flow. That is, we 
have the following. 

Corollary 7 (Persistence). If the vector field r\ is a relabelling symmetry, then the symmetry condition for advected 
quantities C^aft) = persists, provided it holds initially. 

Proof. If the left side of equation (JTTJ) vanishes initially at t = 0, then it continues to vanish for all time t > 0. □ 

Definition 8 (Locally conserved quantities). A locally conserved quantity c(t) follows from the equations of motion 
and satisfies a local conservation law, 

(d t +C uit) )c(t) = 0, (18) 

which has the same form as an advection law. 

Remark 9 (Local conservation laws). One distinguishes between advected quantities and locally conserved quanti- 
ties. Namely, equations for advected quantities are obtained from the action G x V — > V and are independent of the 
fluid velocity. In contrast, local conservation laws involve the fluid velocity because they arise from the equations 
of motion. 

Corollary 10 (Iterated conserved quantities). If the quantity c(t) satisfies a local conservation law H8\) as a result 
of the EP equations of motion for all time, then C v c(t) is also locally conserved for any relabelling symmetry. 

Proof. This follows from replacing a(t) by c(t) in equation (fTT|) . □ 

Remark 11. Iterating this process further is possible, but once a conserved quantity can be expressed in terms of 
advected quantities, iteration does not lead to new information. 

In the following section we adopt the strategy of constructing the symmetry vector fields n in Noether's theorem 
using the advected quantities that they preserve. This is accomplished in the examples below for several representa- 
tive fluid flows in three dimensions. Occasionally, when enough freedom remains in the infinitesimal Lie symmetry 
?7, the local conservation law that emerges from Noether's theorem may be re-substituted into the weak form of 
Noether's theorem to compute an additional integral conservation law. 
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3 Examples 

3.1 Advected density: Conservation of vorticity and helicity 

For the specific case that the mass density a = pdV is advected and other advected quantities are absent, e.g., in 
barotropic fluid dynamics, the symmetry condition (|13|) is 



C v {pdV) =d(v-ipdV) = 0. 
Therefore, by Poincare's Lemma, one may write locally that 

ri-\pdV = d(¥- dec) = curl*- dS, (19) 

for some vector function *. 

For non-trivial topology (on a spherical annulus, for example), we may choose a simply connected patch bounded 
by a simple closed curve C{t) that is transported by the fluid velocity u. We then restrict r\ at each time to the Lie 
algebra of vector fields that leave C(t) invariant. (These vector fields n are tangent to the curve C(t).) This choice 
allows us to define * on the patch enclosed by C(t) for each relabelling symmetry rj. 

Equation (|f 9p the vector field 77 for each relabelling symmetry may be expressed in terms of a vector function 
*, as 

1] = p _1 curl* • V. 

All such vector fields satisfy the advection condition (|T5|) for the density pdV, since 

C v {pdV) = div(p / 9~ 1 curl*) dV = div(curl*) dV = 0. 

We substitute this solution for rj into Noether's theorem and use Corollary [7] (persistence of the symmetry 
relation) to find, 

d /SI 
° = di^ 

dt J v du 

= 77 / -7--dxA?7_l(pdF 
dt Jd p du 

By (HI = — / -— -da;Ad(*-da;) 
dt J v p du 

f d (\ SI \ 1/T , , f 1 51 d J/T J , 

= / tt; -7- -da; Ad(*-dx)+ / -— • dx A — d(* • dx) 
J v dt\pSu J JvpSu dt 

By(H3) = / ^ f-^-.dx') Ad(*-da;)+ / • dx A (- d£ tt (* • dx)) 

l d G^- da; ) +/: " d G£- da; )) A( *- da;) 
l +£ ^) d G^- da; ) A(v,, - da;) - 

Since * is arbitrary, vanishing of the final line gives the weak form of the local conservation of the vorticity 
2-form, 

i^)^'")-'' (20) 

where vorticity is defined as the curl of the specific momentum (momentum per unit mass) . The specific momentum 
is equal to the velocity for Eulcr fluids; so its curl in that case is the usual Eulcr fluid vorticity. 

Remark 12 (Two velocity vectors). Two velocity vectors appear in the computations above: These are the fluid 
velocity vector u in the Lie derivative C u and the specific momentum covector p~ 1 5l/5u in the 1-form p~ 1 Sl / Su-dx. 
These two velocities are the basic ingredients for performing modelling and analysis in any ideal fluid problem. They 
appear together and have separate meanings in the Euler-Poincare equation and throughout the present paper, as 
illustrated in the examples below. 
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Example 13 (Incompressible Euler equations). The incompressible Euler equations have reduced Lagrangian 
where p is a Lagrange multiplier enforcing incompressibility p — 1 . The variational derivatives are 



81 81 \u\ 2 

In this case, the conserved vorticity is 



pu ' T P = — - p - 



curl [ — —] = curlu. 
\p8uj 

Example 14 (Incompressible Euler-alpha equations). The incompressible Euler-alpha equations with p = 1 have 
reduced Lagrangian 

l(u,p) = |(M 2 + a 2 \Vu\ 2 ) +p(l -p)dV. 



The variational derivatives are 



SI •>„ ^ 31 l/i .2 2m |2\ 

— = pu - a 2 V ■ pVu, — = - ( u Y + a 2 Vit 2 - P- 
ou dp 2 



In this case, the conserved vorticity is 



curl — — = curl it V ■ pvu, 

\p8uj p 



which becomes curl(it — a V u) since p = 1. 

In the preserved 2-form d(\& -da;) introduced in equation ([T^|. the vector function *f? is determined (locally) for 
each choice of symmetry vector field r\. Likewise, we have seen that each choice of the vector function \& corresponds 
to a certain relabelling symmetry r\. Consequently, Corollary [7] for the persistence of symmetry and the definition 
of a local conservation law in equation (fT8")) would allow us to replace d('i r • da;) with another conserved 2-form. In 
particular, we may choose the conserved vorticity 2-form and set, 

( 1 81 

d(* • dx) = d I - — • da; 



p 8u 

After this identification, we may draw the conclusion from persistence in Corollary [7] that 

51 



= d^W 77 

d f 1 81 
= — I - — ■ dx Arj -1 pdV 

dt J v pbu 

A. f i *L da;Ad (- — dx 
dt J v p Su \p Su 

dt fx, p Su i^p 8u} 

Thus, the weak form of the local conservation law for vorticity yields a conservation law for the helicity integral, 

f 1 81 (\ 81 \ 
H '= / ZTT.'^ -I- )*V. (21) 



p Su \p8u 

The vector field for the symmetry associated with helicity conservation is 
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The characteristic paths of the vector field tjh may be regarded as vortex lines, and these satisfy the symmetry 
condition ([T2"j) . as a result of the EP equation ©. That is, the characteristic paths of r/n are frozen into the flow of 
the fluid velocity. This means that shifts along these paths are relabelling symmetries and the corresponding Noether 
conservation law is the helicity H in equation pip . In particular, the symmetry associated with conservation of 
helicity is a relabelling of the frozen-in vortex lines. 

As mentioned earlier, the momentum per unit mass is equal to the velocity u for Euler fluids, so that conservation 
of the helicity for Euler fluids may be expressed as 



dt 



u ■ curl u d V = . 



The spatial integral H defining the fluid helicity in ([2"Tj) measures the knottedness, or number of linkages, of the 
vortex lines, that is, lines of ui = curl(p _1 i5Z/(Sit). This fluid helicity indicates the topological complexity of the 
winding of the vortex lines in uj amongst themselves in the spatial domain [3]. Physically, helicity conservation 
arises because the vortex lines are frozen into the flow of the diffeomorphisms and, thus, they cannot unknot. 

Remark 15 (Ertcl's theorem in hydrodynamic notation). We identify the evolution operator in Theorem [S] as the 
familiar Lagrangian time derivative D/Dt, 

dt + £ u (t) = ^ , 

and we express the vector field r/H in equation (|22p in terms of a generalised vorticity vector oj, defined as 

-1 , ( l Sl \ „ -1 r, -i ,/HI 

Vh = P curl — — • V =: p to ■ V with w := curl — — 
\pouJ \pdu 

Introducing this familiar hydrodynamic notation allows one to write the symmetry relation (|16j) in the case for the 
action of the Lie derivative on a function a(t) as 

^(p- 1 ^ • V)o(f) = (p-'ui ■ V)^a(t) , (23) 

which is the usual form of the classical Ertel theorem [10] . For a scalar advected function, a G A , Corollary [7] 
(persistence) yields yet another scalar conservation law, for q = {p~ l u) ■ V)et. 

3.2 Advected density and tracer: Conservation of potential vorticity 

Proposition 16. For the case of two advected quantities a\ = pdV € A 3 , 02 — s G A , the simultaneous solution 
of C v p d V — and C v s = is 

n_lpdV = d(0ds), (24) 

for general <f> G A . 

The proof of this proposition is simple, because a\ is a top form and 02 is a bottom form. 
Proof. Symmetry requires that these two advected quantities satisfy 

77 _l p d V = d(\& • dx) and n _l ds = 77 • Vs = . 



Thus, 

and, hence, 



dsAn_ipdV = (Vs-n)pdV = 0, 
= (Vs • rj) p d V = ds A (77 _l p d V) = ds A d(* • d a;) = ds A d (0 d s) 



□ 



The two advected quantities d(^f ■ dx) and d (<f>ds) both equal (n _l pdV), so they satisfy the same evolution 
equation. In particular, the following advection equation holds 

-d(^d S ) = -d£„(^ds). 
ot 
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Now one may substitute r\ _l pdV = d(</>ds) into the Noether theorem calculation as above and recompute, finding 
this time that: 

d /SI 

= d7\^ 

By® = J- t J v (- p ^- dx ) Ad(0d s ) 

By® = / ^ f-^-d^ Ad(0ds)+ /" -^--d^AC-d^^ds)) 
7x> <« \pou J J v pdu 

^-d (-^- -dx] Ads + £„d [ ■ dx ) Ads, <f>) 
at \pou J \pou J I 

Bys = a 2 = -(|(d(i^. d:E )Ad S ) + /:„(dQ£.dx)Ad S 

£u(t) | | df --^ - dx ) Ads 



9i / V \P Su 

As before, all boundary terms have been dropped in spatial integrations by parts. Since <j> is arbitrary, the final line 
of the calculation above gives the weak form of the conservation law for Ertel potential vorticity (PV) density, 
defined as flQ] , 

qpdV := d ( -dx ) Ads = curl( ) -VsdV. (25) 

\pou J \pouJ 

The corresponding local conservation law is 

^+£ u }(qpdV)=0. (26) 



/ $(q)pdV = 0, (27) 



The arbitrary function (/> in the weak form of the local conservation law for potential vorticity also yields the integral 
conservation law, 

d 

di 

in which $ is an arbitrary function and we used dt(&(q)p) = — V • (&(q)p). 

The vector field for the symmetry associated with PV conservation in (|26l) may be computed from equation 

as 

Vpv = P^ 1 (V</> x Vs) • V, (28) 

and it represents shifts along level sets of s. In particular, the characteristic paths of the vector field r\pv he 
along the level sets of s, which in turn are frozen into the flow of the fluid velocity. This means that shifts along 
the characteristic paths of r\py are relabelling symmetries and the corresponding Noether conservation law is the 
advection of the potential vorticity q in equation (1251) . 

Example 17 (Rotating Euler-Boussinesq equations). The reduced Lagrangian for the rotating Euler-Boussinesq 
equations is 



I = 



/ P~i~ + pu- R- zb + P (l- P )dV, 
Jv * 



where b is the buoyancy satisfying ^| + C u b = 0, and where R satisfies curli? = 2ft. Consequently, the conserve 
potential vorticity may be computed, as follows: 

curl ( -— | = curlfu + R) with curlB = 2n =^> q = (curlu + 2ft) ■ V6. 
\p SuJ 
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3.3 Advected density and flux (2-form): Conservation of cross helicity 

Proposition 18. For the case that a\ = pdV G A 3 and 02 = B ■ dS = d(A ■ da;) G A 2 , the only simultaneous 
solution of C n pdV = and C^B ■ dS = is 

rj-ipdV = B ■ dS. (29) 

Proof. 

Recall n _l pdV = d(* • dx) and identify d(* ■ d x) = B ■ d S . 

□ 

In this case, Noether's theorem implies the conserved quantity 

d / SI \ 

By = 4- f --r- ' dx A B -dS 

at Ju p 011 

= ±j (b. 1 -^L 

dt J v \ pSu 

This is the cross helicity, which is known to be conserved, in particular, for ideal magnetohydrodynamics (MHD) 
[3]. The vector field for the symmetry associated with conservation of cross helicity is, 

Vch = P~ X B ■ V, 

which represents a field of shifts along magnetic field lines. The characteristic paths of the vector field r\cn arc 
magnetic field lines that satisfy the symmetry condition (fT2]) , as a result of the advection equation (|13| for magnetic 
flux B ■ dS. That is, the characteristic paths of t\ch are frozen into the flow of the fluid velocity. This means that 
shifts along these characteristic paths are relabelling symmetries and the corresponding Noethcr conservation law 
is the cross helicity. 

4 Other conservation laws for ideal fluids 

In this section, we first point out that not all fluid conservation laws follow from Noether's theorem, as formulated 
above, by considering the counterexample of magnetic helicity for MHD. We then make a connection between the 
Noether's theorem discussed in this paper, and the Kelvin- Noether circulation theorem discussed in [To] . 

4.1 Magnetic helicity 

The distinction between advected quantities and locally conserved quantities comes back into play, when one 
considers compound advected quantities that arc conserved independently of the motion equation. For example 
advection of the scalar s and the exact 2-form B ■ d S — d(A • d a;) lead immediately to advection of the compound 
quantities, 

ds AB -dS = div(sB)dV and A • d x A B ■ d S = A ■ B d V . 

The former is trivial, because its integral over space vanishes identically. However, the latter is the famous mag- 
netic helicity, whose spatial integral measures the knottedness, or number of linkages, of the magnetic field lines. 
That is, the magnetic helicity indicates the topological complexity of the winding of the magnetic field lines amongst 
themselves in the spatial domain. The preservation of this magnetic winding number is a fascinating property of 
ideal MHD flows [3J, but it does not arise from a Noether symmetry. It arises here as a compound Lagrangian 
quantity whose Eulerian interpretation is deep and interesting. It is beyond our present scope for further study, 
except to provide a counter example to the conjecture that a converse of Noether's theorem might exist for Euler- 
Poincare ideal fluid theories. 
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4.2 Modified vorticity, potential vorticity, helicity and the Kelvin-Noether theorem 

More general conservation laws can be obtained by expanding the set of variables, so that the time variation of the 
quantity a is enforced by a Lagrange multiplier b (known as a Clebsch variable) instead of constraining the variation 
5a. As described in 0EJQTJ, the same Euler-Poincare equations are obtained this way. In this case, Hamilton's 
principle becomes 

/ / Q 



Sciebsch=l l(u,a) + (b,[ — +C u )a)dt, (30) 



with Lagrange multiplier b to be determined. Then, Hamilton's principle yields after a short calculation, 
= SSciebsch = (Jr^ -boa, SuJ + (^5b, ^ +£, 

l-^ t+ C u )b,Sa)dt + [(bM)}l. 

We now consider symmetries of the form 

5u = r) + [it, 77], 5a = 0, 8b = 0, 

for a general relabelling vector field 77 that satisfies 77 + [u, 77] = 0, but is not constrained to be a symmetry of the 
quantity a. Noether's theorem then leads to 

d / 1 61 \ 

Next, we define the (o) operation in terms of the diamond operation by 

boa =: (~ da ^j ®P&V- (31) 

The ( o) operation allows one to express a 1-form density as the product of a 1-form and the advected mass density. 

After a calculation similar to that leading to the result (|20l) , one may write the vanishing of the 77-coefficient in 
the previous variational equation for SSciebsch = as 

57 fj I - Tt { r ' + £ G e) 50 - L G ls ' + c - )b ) 50 <52) 

which is found after substituting ( + Cu) a = 0, as imposed by the 5b- variation. In the loop integrals, the closed 
circuit 7 t moves with the flow of the fluid velocity vector field u. 

Now we have two choices. Namely, we may either eliminate Lagrange multiplier b by using the variational 
equation for b, 

l + £ "> = ^ < 33 > 



or we may keep b as an additional dynamical variable satisfying equation (|33[) . The first choice yields the Kelvin- 
Noether theorem of [T5], and the second choice yields an advection equation for a quasi- vorticity vector field, plus 
the additional equation (|3"3"|) for b. Specifically, in the first choice, the second and fourth terms in equation (|3"2"j) 
cancel, leaving 

d r 1 51 r fl 51' 



dt J yt p 5u ' J lt \p 5a) ^' ^ ^ 

which is the Kelvin-Noether theorem [15] for circulation. 

In the second choice, the third and fourth terms in equation (|32[) cancel instead, thereby leaving the following 
circulation conservation law, 

d / fl 51 b _ \ n 



d t J lt \ p 5u 

or, equivalcntly in vector notation, 



— (f u ■ d x = , with u ■ d x := ( oa) . (36) 

dt J \p 5u p 
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The price for this circulation conservation law is that the Lagrange multiplier b remains and satisfies equation ([53]) , 
instead of being eliminated. However, keeping b as a dynamical variable also has the added value that doing so 
yields a quasi-vorticity Cj involving b that satisfies the advection law for a 2-form under the flow of the velocity 
vector field, u, 

— + C U ) (5'dS) = 0, with £j • dS = (curl u) • d S := d ( — - on . (37) 

ot J \p ou p J 

Moreover, the vector field p~ 1 uj, which is derived via the relation 

p^u-lpdV = ib ■ dS , 

also satisfies the invariance equation (|12[) . namely, 



9 I 



U , p Ul 



= , which means 



ph. + ^ u > 



= 0, (38) 



as demonstrated in equation ([T5|) in the proof of Theorem 2) 

Consequently, keeping the Lagrange multiplier b as a dynamical variable produces an Ertel theorem of the form 
(|23p and yields conservation laws for the corresponding potential vorticity and helicity. 

Moreover, these equations apply for essentially any fluid theory; so keeping the Lagrange multiplier b instead of 
eliminating it affords a certain universality to the formulation. The corresponding conserved potential quasi-vorticity 
q and quasi-hclicity H arc defined by 

q := p^ 1 uj ■ X/ s with H := I u ■ dx A d(u ■ dx) = / <D • curl -1 uj d V , (39) 

Jv Jt> 

in terms of the quasi-vorticity £> in (|37[) and writing g for an advected scalar function, a\ = s, and mass density, 
0-2 = pd V. 

Example 19. As an example, consider the particular case a\ = s and a2 = pdV, when an analogue of potential 
vorticity exists. In this case, the quasi-vorticity 2-form is given by 

<2> • dS = d ((it + p~ l b Vs) ■ da;) = curl (it + p~ x b Vs) ■ dS, 

and the Lagrange multiplier 6 satisfies 

Thus, at the cost of keeping the 6-cquation as ( + £ u ) b = 81/ 8a, we can extend the Kelvin, Ertel and helicity 
theorems for most fluid theories, but these are Noether conservation laws for the Lagrangian corresponding to the 
extended variable set (it, a, b). 



d \ b 1 SI 

u ■ V 



5 Conclusions 

In this paper we showed how to obtain conserved quantities for ideal fluid models that can be obtained from the 
Eulcr-Poincarc equations with advected quantities. The conserved quantities are obtained via Noether's Theorem 
from relabelling symmetries of the Lagrangian, which are generated by vector fields that satisfy the condition 
r\ t — ad„ r] = 0. Fluid theories usually involve advected quantities that evolve according to the equation a t + C u a = 0: 
an advected density is almost always present, and other possibilities include advected scalars such as buoyancy, or 
advected 2-forms such as magnetic flux. In order to have a symmetry of the Lagrangian it is also necessary to 
satisfy equation (fT3|) so that the vector field 77 generates a symmetry of the advected quantity a(t). In Corollary 
[7] we showed that if (fT"3|) is satisfied initially then it is satisfied for all subsequent times. In general, defining a 
parametcrisation of the null space of C v for 77 then leads to evolution equations defined on the dual of the space 
of parameterising functions. Note that this approach is different to that of Section 4 of [15], which does not use 
Noether's Theorem and instead performs computations on the Euler-Poincare equations directly. 

In this paper we considered fluid theories with advected density, leading to conservation of vorticity, as well 
as advected density plus advected tracers, leading to conservation of potential vorticity, and advected density plus 
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advected 1-forms (such as magnetic flux), leading to global conservation of cross- helicity. There are many other 
advected quantities that could be considered, notably advected tensor fields which are used in the theory of ideal 
complex fluids |13j . However, since most fluid theories include an advected density (even incompressible flow, for 
which the pressure is a Lagrange multiplier that enforces that the density remain constant), it is not always possible 
to simultaneously solve all of the constraints on 77 arising from the requirement that 77 generates symmetries of 
all of the advected quantities involved. For example, when a density and a magnetic flux are both present, there 
is only one symmetry and hence only one globally conserved quantity. An interesting class of problems in which 
density is not necessarily present arise in computational anatomy [17j . Here the aim is to find the solution of the 
EPDiff equation (for which the reduced Lagrangian is a functional of u only) which transports one configuration 
of an advected quantity a to another. The advected quantity might be a scalar (for greyscale images), a singular 
measure (for curves and surfaces), or even a tensor field (for diffusion tensor images). Solutions of the EPDiff 
equation are geodesies on the diffeomorphism group; for these solutions to drop to geodesies on the shape space 
corresponding to the chosen advected quantities, all of the conserved Noether quantities must vanish QT|. For 
example, for greyscale images described by advected functions I(x, t), the momentum is constrained to be normal 
to the image: 8l/Su ■ VI = 0. Hence, the Noether quantities define the geometry of the shape space. 

One of the "holy grails" in the field of variational numerical methods is to find Eulerian discretisations of fluid 
dynamics that arise from a variational principle. Amongst other things, this would provide the possibility of discrete 
forms of the Nocther's Theorem described in this paper. One direction that we have previously explored is to try to 
find a discretisation of the diffeomorphism group and to obtain some form of reduction by symmetry; this approach 
has been developed in some detail, making extensive use of discrete exterior calculus, for the case of incompressible 
flows in [55] . In [7] , it was shown that the spatial discretisation of the Lie bracket must satisfy the closure property 
if a reduction is to be obtained; it was also shown how space-time discretisations could be obtained in this case. 
Another direction that we have explored is using Clebsch constraints to enforce the evolution of the back-to-labels 
map [8] . This leads to a multisymplectic formulation of fluid dynamics that can be discretised by a standard recipe 
but reduction (elimination of the back-to-labels map) is not possible after discretisation due to symmetry breaking. 

On the Hamiltonian side, the conserved quantities associated via Noether's theorem with relabelling symmetry 
comprise the Casimir functions. The variational derivatives of the Casimir functions are null eigenvectors of the Lie- 
Poisson Hamiltonian structure that arises from the Euler-Poincare framework upon Legcndre transforming. This is 
explained further in [T6j [TSJ [1] . See also [27l [28] for related discussions and additional references. The conservation 
laws that follow from Noether's theorem for relabelling symmetry of the Eulerian fluid variables generate steady 
flows when substituted into the augmented Lie-Poisson brackets that include the particle labels as functions of time 
and spatial coordinate [12l Q] . 

Thus, Lie-Poisson brackets with the Casimir functions leave the Eulerian fluid variables invariant, but they shift 
the fluid particle labels along steady flows. This fact led to a strategy for proving nonlinear stability of equilibrium 
flows that was first recognized in [2] for ideal incompressible planar flows and was later developed and applied more 
widely in plasma physics in |16j . Likewise, the Eulerian vector fields for relabelling symmetries found here could 
just as well have been obtained by solving for the null eigenvectors of the Lie-Poisson bracket on the Hamiltonian 
side. Critical points of the sum of Hamiltonian and the Casimirs are steady equilibrium flows. The stability of these 
equilibrium flows may be studied by taking a second variation and determining the conditions on the equilibrium 
that would make the corresponding linearly conserved second variation sign definite |16] . 

Remark 20 (Noether's other theorem). As mentioned earlier, Noether's original paper actually contains two 
theorems. The second one is generally regarded as the more subtle of the two. Noether's second theorem leads in 
principle to dependence among the Euler- Lagrange equations (Bianchi identities). However, for ideal fluids, we 
have not found any strictly Eulerian conservation laws in addition those already discussed here by using Noether's 
second theorem in the Euler-Poincare context. 
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